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^ ■ Abstract 

, The Bel-Robinson tensor Ba/s^i, gives a positive definite gravitational en- 

' ergy in the quasilocal small sphere limit approximation. However, there is an 

alternative tensor Vap^u that was proposed recently that offers the same posi- 
CJ ! tivity as Bap^u does. We have found that Vap^u is the unique alternative tensor 

with Bafi^u which implies that these two tensors are a basis for expressions that 
have the desirable non- negative gravitational energy in the small sphere limit. 
In other words, the 'energy-momentum' density according to Bap^y and Va/s^u 
are on equal footing at the same limit. 
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^ 1 Introduction 

in 

^ ! The Bel-Robinson tensor B^p^j^v has many nice properties. It is completely symmetric, 

O \ completely trace free and completely divergence free. It is usually regarded as being 

related to gravitational energy. In particular, the gravitational energy-momentum 
density in the small sphere vacuum limit is generally expected to be proportional to 
. the Bel- Robinson tensor. This expectation is related to the requirement of energy 

\ positivity pj. 

_c3 \ It should be emphasized that, generally speaking, a positive energy proof for a 

quasilocal expression is not easy. Here we consider specifically pseudotensor expres- 
sions. In fact, the quasilocal methods are not fundamentally different from pseudoten- 
sor methods [2j. The gravitational energy expression in a small region limit can be 
investigated through the pseudotensors. Normally, the expansion of a pseudotensor 
expression up to second order can be represented by certain tensors Bap^u, Sap^u and 
Kapfiu [31 H] . In other words, quasilocal expression is a fancy name for a pseudotensor. 
Even though a pseudotensor is not a tensorial object, this does not imply that it is 
useless. The second order expansion expression provides guidance as to whether the 
gravitational energy expression is positive or not. More precisely, a negative quasilo- 
cal gravitational energy expression on a small scale definitely guarantees that it can 
be negative on the large scale. Conversely, a positive energy expression in the small 
region implies that there may be a chance to obtain positivity in a large region. 

However, it is natural to question whether studying this kind of quasilocal formu- 
lation has physical significance. It is well known that the gravitational energy density 
cannot be detected at a point because of the equivalence principle (see section 20.4 
in Fortunately, the quasilocal idea has physical meaning, i.e., the gravitational 



energy density is well defined at the quasilocal level theoretically [SlEllTj. Practically, 
pseudotensors can be used to calculate the tidal heating |H] (e.g., Jupiter and lo) as 
well as using the quasilocal formalism [9]. Moreover, from the last decade, there are 
many researchers who believe that finding a good quasilocal expression (especially 
one that is locally positive) is meaningful and worthwhile, see e.g., [H [10] and the 
many references contained therein. 

In the past, the Bel- Robinson tensor has been considered to be the only tensor 
which contributes positivity in the small sphere limit. However, we recently found 
another tensor [3], Vap/^u (defined in (!29|) below), which is also quadratic in the cur- 
vature, and which enjoys the same positivity properties as Bap^y. More precisely, 
the associated 4 momentum in a small sphere is a Lorentz-covariant future point- 
ing non-spacelike vector (see 4.2.2 of pQ). For short we call this property casuality. 
Furthermore, we found that Bap^u and Va/3fj.u are a basis for expressions which have 
the desirable non- negative gravitational energy in the small sphere vacuum limit. As 
we found, Va/3fj.u fulfills the weak energy condition which gives positive energy in the 
small sphere limit, this is then sufficient to argue that the 'energy- momentum' density 
according to Bap^u and Vapf^u are on an equal footing in the same limit. We will prove 
in section 3 that Vai3f_iu is the unique alternative tensor that has the desirable positive 
energy property in the small region limit. 

In this work, we examine some properties of Va(5fj.u and some other quadratic in 
curvature tensors, Sap^u, K^p^u and Wapf^u, which have shown up in the expansion 
of energy in the small sphere limit. For instance, it has already been shown that 
Va/Sfiu does not have the dominant energy condition [3] . We also found another tensor 
V^p^^, which is not restricted to the pseudotensor conservation of energy- momentum 
requirement, but does satisfy the weak energy condition requirement. 

Here we are concerned with finding a suitable form for a pseudotensor (equiva- 
lently, a quasilocal Hamiltonian expression, see [TT]) in the small region limit. For 
the zeroth order term, the pseudotensor gives the mass density as the equivalence 
principle demands. Moreover, we also need to consider the ADM mass at the spatial 
infinity. Combining these two constraints [3] can confine a suitable pseudotensor ex- 
pression. Finally, the non-vanishing second order terms contribute the gravitational 
energy- momentum in a small region limit; these terms are quadratic in the curvature 
tensor. 



2 Quadratic curvature tensors 

There are three basic tensors that commonly occur in the gravitational pseudotensor 
expression [H [T2] 
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where = Rpr^^R''^^'^ ■ Some properties of S'^/j^j, and Ka/s^u [3] that are eas- 
ily verified using the above definitions and the well known vacuum identity [13] 

3 

Sa|3^J.u = S(^al3)(pu) = 'S'(^!/)(a;9) 5 Saji^^ = -(^q^R , Saf^p^ = 0, (4) 

3 

Note that unlike Bajs^u, both SaiSfiu and K^js^y are neither totally symmetric nor 
totally trace free. Obviously, (jl]) and ([5]) already indicate that Sap^iv and Kap^u 
do not have the completely trace free property. For the non-completely symmetric 
property, one can verify this by using Petrov type D [14j. In particular, we found a 
case where 

'S'oOll 7^ 'S'oiOl, -f^OOll 7^ -^'^OlOl; (6) 

where we have evaluated them using the electric part Eat and magnetic part Hab, 
defined in terms of the Weyl tensor [15] as follows: 



Eab '■— CaObO, Hab '■— *C*aOfeO- (7) 

In order to appreciate the nice properties of the Bel- Robinson tensor, we compare 
some components of SapfMu and Kapf^u. For the analog of the electromagnetic stress 
tensor energy density, 

-DOOOO — ^ab^ + nabn , {6} 

Smm = '^{EabE"-^ — HabH"'^), (9) 
i^oooo = -EabE'^' + ^HabH'^'. (10) 

Likewise for the momentum density (i.e., the Poynting vector) 

Boooi = '^'^ijkE^'^H^di 5*0002 = 0, i^oooi = '^eijkE^'^H^d- (H) 
Finally, the stress, 

Booij = 5ij{EabE°''' + HabH"'') — 2{EidE/ + HidH/), (12) 

2 \6,,{-EabE''' + HabH'^') + 4{-E,dE/ + HuH/)] , (13) 



Si 



OOij — 

Koo^J = 6,j{EabE'''-3HabH''')-AE,dE/. (14) 

From the above comparison, it is clear that the Bel-Robinson tensor indeed has the 
best analog with the electromagnetic stress tensor T'^". In detail, in Minkowski coor- 
dinates (r, X, y, z) the components of the electrodynamic stress tensor are 

T°° = ^ (E.E" + 5,5«) , (15) 

= d'^ejabE^B^ = {Ex B)\ (16) 
1 
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S^i {EaE"" + BaB") - 2 (e'E^ + B'B^)] . (17) 



where E and B refer to the electric and magnetic field density. 

Using a Taylor series expansion, the metric tensor can be written as 

9.p{x) = QaM + d^g^p{Q)x^ + \d%g^p{Q)x^x^ + .... (18) 
At the origin in Riemann normal coordinates (RNC) 

^a/3(0) = d^iQapi^) = 0, (19) 

~3(9^j^(?a^(0) = Rafi/Su + RauPfj,, —"idyV^ apiS^) = R^ ajSv + R^ I3av ■ (20) 

For the quadratic curvature tensors, there are 4 independent basis [IE] expressions, 
we may use 

BaP^u • RaXfiaRp u ~l~ RaXuaRp fi Bo^p^i, + -g^pg^,Yi\ (21) 

o 

SaPfiu • RafiXcjRpu RauXa Rpfi — SaPfiu — -^gapgfiu^ , (22) 

A A 3 2 

Kap^u '■= RaXPuRfj. + RaXPcjRv fi'^ = + -^gapgfiu^ , (23) 

O 

Tapfiu '■= —-gapgfj.u^^ ■ (24) 

O 

These four tensors are manifestly symmetric in the last two indices which means 
Map^,y = MapiiMu)- Then it automatically impjy M^p^^^) = M(^ap){i,u)- Moreover, 
it also naturally turns out that M(^ap){iJ.u) = ^(^i/)(a/3)- Explicitly, they fulfill the 
symmetry Map^u = M(^ap){fiu) = M(^^^)(_ap)- Although there exists some other tensors 
different from Bap^u, Sapfiu, Kapfj.u and Tapfj,u, they are just linear combinations of 
these four. For instance 

TafiPu + T^njjp^ = B^Pnu + ~Sap^i, — K^p^u + '^T^p^u- (25) 

The above identity can be obtained by making use of the completely symmetric 
property of the Bel- Robinson tensor. Using ( !25|) . we can rewrite the Bel- Robinson 
tensor in a different representation [16j : 

1 5 1 

BaPnv = —-SaPiMU + KaPfj,u + -gapg^iv^ — -{gaij,gpv + gaudPti)^ ■ (26) 

Z o o 

This equation will be used in section 3. 



3 An unique alternative non-negative gravitational 
energy tensor in small sphere limit 

3.1 Proof the unique alternative non- negative energy tensor 

Using RNC Taylor series expansion around any point (i.e., at any preselected point 
we consider a small coordinate sphere in RNC, see e.g., PQ EJ [IZl [IS]), consider all the 



possible combinations of the small region energy-momentum density in vacuum. In 
the neighbourhood of any preselected point in RNC, the pseudotensor then has the 
form [T9l 



3^ 

5 



(27) 

where k = 8ttG/c^ (here we take units such that c = 1 for simplicity) and ai to 
are real numbers. Here Ga^ is the Einstein tensor, but we will consider the vacuum 
case, so 

GJ = kTJ = Q. (28) 

Then the first order linear in Ricci terms (9(Ricci, x) vanish. The lowest order non- 
vanishing term is of second order, compared to this in the small sphere limit we can 
ignore the third order terms 0{x^). From now on, the second order term will be kept 
but the others are dropped. The essential purpose of the present paper is to prove 
that Bapf,„ and Va^f.^, where [3] 

ValSfiu '■= SaPfjiu + Kai3^u, (29) 

are a basis for positive gravitational energy in the small sphere limit. There are two 
physical conditions which can constrain the unlimited combinations between Bap^ui 
Sa^^vi Ka/SfMu and T^p^u- The first one is the conservation of the energy-momentum 
density and the second is the positive gravitational energy in the small sphere limit, 
(i) First condition: energy-momentum conservation. Consider fl27|) as follows 

= dpt^^ 



2 [aiB 



= -(ai — 2a2 + 3a3 — a4^)gai3X^^'^ ■ (30) 
Therefore, the constraint for the conservation of the energy-momentum density is 

04 = ai — 2a2 + 303. (31) 

Although there are an infinite number of combinations which can fulfill the above 
constraint, it has removed one degree of freedom. As each single tensor of Bap^ui 
SaPnv^ Kai3^u OT Tap^y cauuot Satisfy the conservation requirement, but a linear com- 
bination of them can. One can simplify the situation by eliminating Tap^y which is 
absorbed by Bai3nu, Sap^iu and Kap^y. Then there are only 3 basis tensors left. Thus 
one can rewrite f l27|) as 

2k = aii^BoP ^y + TfJ^ + (X2{SqP ^y — 2Ta^ -\- dzi^Kj^ + r^Tj^ x^x 



= {aiBj^y + a2Sj^y + a^Kj^y) x^x\ (32) 

Paying attention to [HI [20], when we consider all the possible expressions for the 
pseudotensors which including the fiat metric, there does appear a linear combina- 
tion of these three tensors. We defined Ka^Lv [12] just for convenience and without 



any physical reason. In the beginning, it seems an interesting and even a mysterious 
combination that why only Bap^u, Sap^u and K^js^u always showed up in the expres- 
sion. Nothing more, nothing less. Now, we discovered that it is not an accident but 
becomes a necessary requirement because only Bap^u, Sap^u and Kapf_iu can satisfy 
the condition of the energy-momentum density conservation. 

(ii) Second condition: non-negative gravitational energy in the small sphere limit. 
The purpose of the pseudotensor is for determining the gravitational energy-momentum, 
the associated energy-momentum can be calculated as 



Jt=o Jt=o 



= t\i^— j r^d^x = t\i' (33) 

where l,m = 1,2,3. Using this calculation method, the energy-momentum in the 
small sphere limit for (|32ll becomes 



5 

= i-E, P) = (ai5^o/' + a^S^^i' + a^^K^^i') . (34) 

The 'energy- momentum' values associated with Bap^u, Sap^u and K^js^v are propor- 
tional to 

Bfj,Qi = {EabE""^ + HabH"-^ ,2ecabE°''''H^d)i (35) 
5^0^' = -l^iEabE'^' -HabH-\Q), (36) 
Kfj.01 = B^Qi — 5^0/'- (37) 

Looking back at (1321) . we are interested in the positive gravitational energy within 
a small sphere limit, the Bel- Robinson tensor already satisfies this condition. Precisely 

Bqqi = EabE"'^ + HabH"^^ > 0. (38) 

The rest of the job is to find the coefficients 02 and 03. Using fl37|) . rewrite flM|) as 

,5 

\a, + a3)B^o/ + {a2-as)Sf,o/ • (39) 



P =--L- 
^ 60G 



Equation (136!) shows that 5*^0/' cannot ensure positivity, since we should allow for any 
magnitude of H-Eabll and H-ffafell- In other words, for Sap^^u the sign of the 'energy' 
density is uncertain. Therefore the only possibility for (15^ to guarantee positivity is 
when Oi + 03 > 10|a2 ~ c^sl- However, if we consider that the momentum should be 
future pointing and non-spacelike (i.e., inside the light cone such that —Pq > \P\), 
the unique requirement for (1391) to assure non-negative is when 02 = 03. In other 
words, causality. Moreover, using f l2^ and ([ST]), we obtained 

V^oi=S^o/ + K^oi=B^oi. (40) 

Consequently (!MI) becomes 



5 5 
P^^ = (ai^^oz' + a2V^oi^) = + a2)5^o/- (41) 

Hence the proof is completed. Indeed Bap^u and Vap^u are a basis for expressions 
which have non-negative gravitational 'energy' density in vacuum. 



3.2 Physical meaning of the completely traceless property of 

Bf3^iy and Vap^iy 

For the quasilocal small sphere region, there are four fundamental quadratic Weyl 
curvature tensors. We know that MaPiiu = M(^ai3){fj.u) = for M^p^u ^ 

{Bap^u^ Sap^ivi Kai3^u, Tajs^y}. lu Order to check the completely trace free property of 
the linear combinations of Bap^u, Sa/Sfiu, Kai3^u and Tapfj,^- We only need to consider 
two cases. 

Case (i). Consider the trace on the first and third indices: 

(iiB°'^au + 0-28°' fiau + asK^^au + CL^f'^^au = o (c^i " + 803 — a4,)g^^Il^. (42) 

o 

Case (ii). Consider the trace on the first pair: 

For completely traceless, fB2|) and fH3l) have to vanish at the same time 

ai — 2a2 + Sa-s — = 0, (44) 
01 + 02 — 04 = 0. (45) 

Notice that (jSj) is the same constraint for the energy-momentum conservation ex- 
pressed in fl3T]) . This means that one of the mathematical trace free condition turns 
out to be one of the physical criteria. Solving the above two equations, we recovered 
the same requirement for the gravitational energy-momentum (i.e., casuality) which 
was indicated in fl5^ . explicitly 02 = 03. Moreover, using this totally traceless prop- 
erty at the quasilocal small sphere limit shown in flH|) and fH^ . we recovered the 
same result: that there are two tensors that generate the basis mentioned in section 
3.1. The proof is follows 

(^iBaPfiu + Cl2Sal3fj.u + (^S^afSfiu + O-iTap^u 
= ai{Ba/3f_iu + Ta^u) + 0'2{Sal3nv + K^p^iu + TaPfj,y) 

= CLlBaP^u + a2VaPf,^, (46) 

where 

+ T^p^u. (47) 

The completely traceless property turns out to be a remarkable result because it is not 
just a mathematical property, but it reveals some physical meaning and conditions. 
Namely, the conservation of the energy-momentum and casuality. 

In other words, we have discovered necessary and sufficient conditions. From con- 
sidering the quasilocal small sphere limit, we find the completely traceless for these 
two fourth rank tensors Bai3^u and Va^^u guarantees the fulfillment of the energy- 
momentum conservation and casuality. Conversely, if any fourth rank tensor satisfies 
the conservation of energy-momentum and casuality, then it must be a linear combi- 
nation of Bafj^u and Va^^^u- 



3.3 Counting the independent components of Bap^v and Vap^v 

As Bafj^v and Vap^v are different but share the same gravitational 'energy-momentum' 
density, then one may interested to know how many non- vanishing independent com- 
ponents Bap^y and Vap^v have. Using f l26|) and ( 129|) . here we write the ahernative 
relationship between B^js^v and Vajs^i, [3] 

+ W^p^^, (48) 

where we have defined 

3 5 1 

Basically, these three tensors {B, V, W) are fourth rank, in principle they could have 
256 components. However, after considering the symmetry properties, they only have 
a relatively small amount of independent components. It may be important to do 
this simple counting because it reduces the workload (i.e., computer algebra) when 
one calculates all the components of these three tensors. 

First of all, we count the number of components of Bapf^u- In principle, as Bai3^u 
is completely symmetric, by explicit examination it reduces to 35 components. There 
is a formula that directly gives this number. A th rank totally symmetric tensor in 
n dimensional space has CJ^~^^~^ components. For our case 

= 35. (50) 

Since Bap^u is completely trace free, there are 10 constraints which can be replaced 
by the other components. Finally, we only left 25 independent components for Bai3^u 
(see 0). 

Secondly, we count Va/s^u- Tensor V^^^y does not have the completely symmetric 
property, but fulfills some certain symmetries Vap^u = V(^ai3)(fj.u) = V{iiu)ia/3)- In prin- 
ciple, this reduces Va/s^u to 55 components. However, when we consider the totally 
trace free property of Vai3^u (but not completely symmetric) , then there are two extra 
constraints need to be taken into account 

Va^.u = 0, V^%„. = 0. (51) 

Hence we have 55 — 10 — 10 = 35 independent components for VajSfiu- 

Finally, for the completeness, we count Wapfiv Note that Wajs^iu and are 
similar. In detail 

In principle, there should be at most 35 components. However, one must consider 
the extra constraint 

W^p^, + W^^,p + W^^fi^ = 0. (53) 
Finally we have 35 — 25 = 10 independent components for Wap^y 



3.4 Physical application for Vap^v 

The physical apphcation of is similar to B^p^j, in the small region limit. For 

instance, the Einstein pseudotensor does not have a suitable positivity expressions as 
the sign of the energy density is uncertain. Recall the result for Einstein in vacuum 

m 

2k E J = 4(45/^- - Sj,,)x^'xr (54) 

io 

Referring to (J33l) . the corresponding gravitational energy is 

5 5 

{ABm'-Sm') = -^{yEa,E'''>-?>Ha,H'''>). (55) 



However, the Papapetrou pseudotensor gives a value which is a linear combination of 
BaPfjiv and Vap^v at the second order [3l [H]. This indicates that there is a chance to 
obtain a positive energy expression in the large region. Recall the result for Papa- 
petrou in vacuum: 

2«:P"'^ = i(45"'^^, - V''^^,)x''x\ (56) 

y 

Similarly, the corresponding gravitational energy from ( 133|) is 

5 5 
^0 = -^(4^00/^ - ^00/) = -^{E^.E'^' + H,,H^'). (57) 

Before we go on to study any further, however, there comes a question whether 
VapfMu and Bajs^u are totally equivalent? Although they have some components that 
are exactly the same (e.g., V^ooo = -B^ooo), we find that Va^^i, and Bap^y are indeed 
different tensors [3]. It is easy to clarify that these two are different, since they are 
defined by different fundamental quadratic curvatures, explicitly 

Bafj^u — Bap^v ~\~ Ta^y, (58) 

In particular, Vap/j^u is completely trace free but not completely symmetric [3]. The 
following lists some properties 

Vapij.u = V(^aP){^lu) = V(pu){aP), Vap^^ = Q = Va^p^ , (60) 
VboOO = VoQl = Vm^l = Kni™' = EabE"'^ + HabH""^ = -BoOOO; (61) 

^mooo = = Vfj.io'' = (EabE"-^ + HabH"^^-, 2ecabE"''^H^d) = B^qiK (62) 

It is known that Bap^v has the dominant energy property |2T1 [22] 

Bappiv w" wl > 0, (63) 

where Wi,W2,W3, are any future-pointing causal vectors. While Vap^u only satisfies 
the weak energy condition and Wap^u fulfills none of them. However, we found W^p^u 
has some interesting properties. A simple computation using ([2]) shows that 

W^p^^u'^tHn" = 0, W^p^^u'^u^u^u" = 0, (64) 
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Ell + Hli) 
{Ell — E22) 

_ +{Hll — -f^22)^ _ 
{EI2 + HI2) 
EI2 + -^^22) 





Table 1: Five different Petrov types 



where t is a timelike unit normal vector and u can be timelike or null. Looking at 
( H8ll . Vapfj,y contains more information than Bajs^u, however it seems that Bapfj^u is 
the important part of Vap^v and Wap^y is a kind of gauge freedom (i.e., it has no 
important physical effect). 

A physical reasonable energy-momentum tensor should fulfill certain energy con- 
ditions. In particular, the local energy density measured by the observer with a 
4-velocity should be non-negative. This energy condition must be true for all time- 
like unit normal vectors ^23j. In fact, we found Vq/j^j, has the non-negative 'energy' 
property 

v^p^XtHn^ = B^p^xtHn^ = EabE^"' + h^.h'^' > o. (65) 

Looking at this weak energy condition and from the reason of the continuity [23], 
the above inequalities must still be true if the timelike vector t is replaced by a null 
vector V. Indeed, we found 

Va^^yv'^v^v^v'' = B^p^yv'^v^v^v" > 0. (66) 

Here we take a simple test, with 3 different cases, to evaluate the inequality of (166!) 
and get some idea what is the value can be. Without loss of generality, consider 
3 simple cases from the unit normal null vector vi = /c(l, 1,0,0), V2 = /c(l, 0,1,0), 
vs = (1, 0, 0, 1) and is a constant. We found 

Bc.pf.yvy^iv'^vl = k'^{Boooo + ABoooi + QBqqu + ABom + Baa), (67) 

where i = 1,2, 3. More precisely, we have calculated the results referring to the five 
different distinct Petrov types [H] as shown in Table [H 

Note that different Petrov types corresponding to different values simply because 
they are evaluated from different cannonial frames. Likewise, it is not surprising that 
different frames associated with different values from the same superpotential such as 
the Freud superpotential [T2l [T6[ ITS l [24] (i.e., using holonomic frames or orthonormal 
frames). The statement is correct according to [23] for Vaj^^y, which is based on the 
fact that Baj3fjLv has the dominant energy property. 

Following from (l33l) . the energy-momentum density for Vq/j^j, in the small sphere 
limit is ^ ^ 

2t,P^ = ^ - r°^o°) = -^^omoo. (68) 



Or, more covariantly, 

2k = V^^^Xf^ft^ (69) 

where 

and it should be recalled that 

EabE^^ + HabH-^ > \2e,abE'"'H\\. (71) 

The physical meaning (non-spacelike energy-momentum) is here simpler and clearer 
than that of the dominant energy condition (!63l) . Obviously V^^^y can play the same 
role as -Ba^^ti/, it ensures a causal 4 momentum in the small sphere limit. In other 
words, the 'energy-momentum' density according to Bap^u and are on equal 

footing at the small sphere region limit. Complementary, the Bel- Robinson tensor is 
no longer the only tensor that have the unique preference for achieving the causal 4 
momentum in the quasilocal small region, but Vap^u can play the same role and it 
becomes the unique alternative choice. 

3.5 Positive energy for the general fourth rank tensor V^^^^^ 

Prom the technical point of view, if we are just interested in positive energy and relax 
the restriction on the pseudotensor constraint, which means the conservation of the 
energy-momentum, there are an infinite number of combinations that have the weak 
energy condition, not including Bajs^y. We define 

where 5,^1,^2,^3 are real numbers and ti + t2 + ty, = 1. Obviously, the energy- 
momentum contribution for Sa/3^u can be ignored according to (122]) . Explicitly 

S^p^uu'^tH^t'' = = S^p^^u^^u^uV. (73) 

On the other hand, 

f^p^uu''u^uV = f^^p^u'^u^uV = f^^p^u'^u^uV. (74) 
Once again, here u can be timelike or null. Using flS^ . rewrite fl72]) 

^alSfiu ~ '^aPuv + (-5 — l)^'^^^,^ + {ti — l)Ta/Sfj,,y + t2Ta^/Su + HTavP^i- (75) 

Based on the weak energy condition and continuity property, we found 

= V^p^.u'^u^u^u' = B^p^.u'^uf'u^u^ > 0. (76) 

This illustrates that there does exist an infinite number of combinations which have 
positivity if we exclude the conservation of the energy-momentum requirement ac- 
cording to the pseudotensor restriction. 



Furthermore, in order to obtain the dominant energy condition, the Bel-Robinson 
tensor is the unique tensor that has the suitable combination from the four funda- 
mental quadratic curvature combinations, namely from fl2Tl) to fl2^ . As a matter of 
fact, BapfMu has more nice properties than the other quadratic curvature combinations 
generally (e.g., SapfMu and Kap^u)- In particular, Bapfj^u possesses the completely sym- 
metric property. However, concerning the gravitational energy at the small sphere 
limit, we found that Vap^v is the unique alternative choice to compare with Ba/s^u- 

4 Conclusion 

Using the four fundamental quadratic curvature tensors, we constructed all the pos- 
sible combinations in the quasilocal small sphere region expression. We recovered 
that Bai3fiu gives a definite positive gravitational energy (more previously a causal 
4-momentum) in the small sphere limit approximation. However, we found an unique 
alternative, the recently proposed tensor V^/j^j,, which also contributes the same non- 
negative gravitational energy density at the same region limit. Based on the two 
physical conditions: energy-momentum conservation and casuality. We found that 
these two tensors can be classified as a basis for expressions which have the desirable 
non-negative gravitational energy in the small sphere region. In other words, Bap^v 
and Vaj3ij,u are on equal footing in the small sphere limit. This means that if we obtain 
Bap^iu or Vaj3^v from the gravitational expression at the small scale, either of them 
is good enough to search whether the expression is positive or not at the large scale. 
For example, the Papapetrou pseudotensor can be a good candidate to study the 
positivity energy expression, as it is proportional to the linear combination of Bap^y 
and Vap^v at the second order evaluation. 

We found that only a linear combination of B^js^u and Vap^v satisfies the energy- 
momentum conservation and casuality physical conditions. Remarkably, the com- 
pletely trace free property for both Bap^u and Vap^v demonstrates the same two 
physical requirements. It may be interesting and also from the practical reason to 
count the independent components for Bap^u, Vap^y and Wap^iu- We found that there 
are 25 independent components for B^p^v, 35 components for Vap^u and 10 for Wap^y 
Moreover, we found that W^p^u, associated with Vapf^u, behaves as a kind of gauge 
freedom. Furthermore, relaxing the restriction of the energy-momentum conservation 
requirement for the pseudotensor, V^p^^, demonstrates that there are an infinite num- 
ber of ways to obtain positivity, namely the weak energy condition. For the conserved 
expressions B^p^u satisfies the dominant energy condition while Vap^u does not, but 
does fulfill the weak energy condition. 
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